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Abstract. The scalar version of the pathway model of Mathai (2005) is 
shown to be associated with a large number of probability models used in 
physics. Different families of densities are listed here, which are all connected 
through the pathway parameter a, generating a distributional pathway. The 
idea is to switch from one functional form to another through this parame- 
ter and it is shown that one can proceed from the generalized type-1 beta 
family to generalized type-2 beta family to generalized gamma family. It is 
also shown that the pathway model is available by maximizing a generalized 
measure of entropy, leading to an entropic pathway, covering the particularly 
interesting cases of Tsallis statistics (Tsallis, 1988) and superstatistics (Beck 
and Cohen, 2003). 

1 Introduction 

When deriving or fitting models for data from physical experiments very 
often the practice is to take a member from a parametric family of distri- 
butions. But it is often found that the model requires a distribution with 
a more specific tail than the ones available from the parametric family, or 
a situation of right tail cut-off. The model may reveal that the underlying 
distribution is in between two parametric families of distributions. In order 
to create a distributional pathway for proceeding from one functional form 
to another a pathway parameter a is introduced and a pathway model is cre- 
ated in Mathai (2005). Section 1 defines this pathway model. In Section 2 
it is shown that this model enables one to go from a generalized type-1 beta 
model to a generalized type-2 beta model to a generalized gamma model 
when the variable is restricted to be positive. More families are available 
when the variable is allowed to vary over the real line. Corresponding mul- 
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titudes of families are available when the variable is in the complex domain. 
This is a pragmatic approach to the distributional pathway. Mathai (2005) 
deals mainly with rectangular matrix-variate distributions in the real case 
and the scalar case is the topic of discussion in the current paper because 
the scalar case is connected to many problems in physics. Further, in Sec- 
tion 3, through a generalized measure of entropy, an cntropic pathway can 
be devised that leads to different families of entropic functionals that pro- 
duce families of distributions by applying the maximum entropy principle 
(Mathai and Haubold, 2007). This is a theoretical approach to the distri- 
butional pathway. Section 4 draws conclusions from the development of the 
distributional pathway. 



2 The pathway model in the real scalar case 

For the real scalar case the pathway model for positive random variables 
is represented by the following probability density, for purely mathematical 
statistical reasons, and to unify and extend the compilation of statistical 
distributions in Mathai (1993): 

f(x) = c x^-^ [1 - a(l - a)x^] x>0, (1) 

a > 0, 5 > 0, /3 > 0, 1 — a(l — a)x^ > 0, 7 > where c is the normahzing 
constant and a is the pathway parameter. The normalizing constant in this 
case is the following: 



c = — r -, for a < 1 (2) 



r(i)r(,^ + i) 
r(i)r(^-l) 

S {af5)i 



for - ^ > 0, a > 1 (3) 

a — 1 



for 1. (4) 



Wc note here with respect to (1) that in information theory and mathematical 
statistics, parameter interpretation is less important than in special cases 
under the interpretation of a physical model or principle (Kaniadakis and 
Lissia, 2004). 
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2.1 Generalized type-1 beta family emanating 
from equation (1) 

For a<lor— oo<a<l the model in (1) remains as a generalized type- 
1 beta model in the real case, where the right tail is cut-off. In a series 
of publications, see for example Mathai and Haubold (1988) and Haubold 
and Mathai (2000), we considered various modifications to the Maxwell- 
Boltzmann approach to the nuclear reaction rate theory and considered the 
situations of nonresonant reactions and cases of depleted tail and tail cut- 
off. It may be observed that all those situations are covered by the pathway 
model in (1). Before proceeding with the properties let us look at some 
special cases first. For a — 0, a — 1, S — 1 with P replaced by /3 — 1 we have 
the regular type-1 beta density, namely, 

fi{x) = Cix''-\l-xf-\Q<x<l, (5) 

where ci = r{^)r{l) ■ ^ ~ 1)7 = 1)<^ — — 1 we have TsaUis statistics 
(TsaUis, 1988; Gell-Mann and TsaUis, 2004) for a < 1, namely, 

f2{x) = C2 [I - {I - a)x]T^ , 0<x<{l-a)-\ (6) 

where Ci is the corresponding normalizing constant. For a = 0, a = 1 in 

(1) we have the power function densities for 7 = 1 and P — respectively, 
namely, 

/3(x) = C3X^-\ 0<a;<l (7) 

and 

/4(a;) = C4 (l-xY, 0<a;<l, (8) 

where C3, C4 are the corresponding normalizing constants. Note that Pareto 
densities (Brouers and Sotolongo-Costa, 2004; Shahzi, 2007) also come from 
the power function models. Further, if 7 = 1, 5 = 1, a = 1, a = 0, /3 = in 
(1) we have the uniform density, 

Mx) = 1, < x < 1. (9) 

Before concluding this section a remark on Bose-Einstein statistics (Ijiri and 
Simon, 1975; Aragao-Rego et al. 2003) in physics is in order. Naturally one 
would expect the Bose-Einstein density to belong to the type-1 beta family. 
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If we take the parameters 7 = 1 and /3 = in the density in equation (5) 
then naturally the function goes to 



g{x) = c* 



1 



1 — X 



and obviously g{x) cannot make it a density in < .t < 1 unless x is 
bounded away from 1. This can be achieved by a transformation of the 
type x — exp (— t — wy) for w > 0, e* > 1. In this case g{x) goes to 



9i{y) 



w 



-, < y < 00, w > 0, e* > 1. 



"1 ^t+wy _ Y ' 

This is nothing but the Bose-Einstein density which is given by 



My) = ce- 



1 



-1 + exp {t + wy) 



, w > 0,0 <y < 00, e* > 1. 



(10) 



(11) 



Remark 2.1 Thus Bose-Einstein density is a limiting form of a type-1 beta 
density of equation (5) where the variable x is transformed by the above 
inverse transformation where the normalizing constant cq cannot be obtained 
through gamma functions but it can be evaluated by appealing to partial 
fractions technique and then taking logarithms and it can be easily seen to 
be the following: 



C6 = w 



In 
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2.2 Generalized type-2 beta family emanating 
from equation (1) 

When a > 1 in (1) we may write 1 — a = —{a — 1),q; > 1 so that f{x) 
assumes the form, 

f(x) ^cx^-^[l + a{a-l)x^]-^, x>0 (12) 

which is a generahzed type-2 beta model for real x. Beck and Cohen's su- 
perstatistics (Beck and Cohen, 2003; Beck 2006) belongs to this case (12). 
For Q; = 2,a = l,5 = l,/3 — 7>0we have the regular type-2 beta density, 
namely, 

f7{x) = C7 x^-\l + x)-^, x>0. (13) 
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In fj{x) for X = ^F, m, n = 1, 2, . . . , 7 = f , /? = f + f we have the F 
density or the variance-ratio density, namely, 

fs{F) = C8 Ff-\1 + !^F)-(f F > 0. (14) 

Th 

For 7 = 1, a = 1, 5 = 1, /3 = 1 in (12) we have TsaUis statistics (TsaUis, 1988; 
Gell-Mann and TsaUis, 2004) for a > 1, namely 

Mx) = cg [1 + (a - l)xr^, X > 0. (15) 

For 6 = 2, —00 < X < 00, a = j;,a = 2, P = in (12) we have the Student-t 
density (Gheorghiu and Coppens, 2004), namely 

/io(x) = cio(l + — ) ~, -00 < X < 00. (16) 

In /lo for 1/ = 1 we have the Cauchy model, namely, 

fn{x) = cii{l + x^)~'^, -00 < a; < 00. (17) 

In if 7 = 1, /9 = 1, a; = e^"*"™ for 77 > 0, e 7^ 0, < y < 00 then we have 
the Fermi-Dirac density (Aragao-Rego et al., 2003) 

fM = ci2[l + exp (e + r)y)]-\ < y < 00. (18) 

Rerricirk 2.2 Observe that /12 is a limiting form of the ordinary type-2 beta 
model with /3 — 7 = in /y. In this case the normalizing constant cannot be 
evaluated with the help of gamma functions but C7 can be evaluated through 
partial fractions and then appealing to logarithms. 

In /y if we transform x to y such that x = e^, —00 < y < 00 then we obtain 
the generahzed logistic and related models 

fisiy) = cismi + e^]-^ -00 < y < 00 (19) 

which are applicable in many areas of statistical analysis, see for example 
Mathai and Provost (2006). Before concluding this section one more remark 
is in order. 

Remcirk 2.3 For the generahzed type-2 beta model, that is the model in 
(1) or the model in (12) for a > 1, x and ^ belong to the same family of 
distributions. Hence we could have replaced x^ by x~^ in the densities in 
equations (13) to (19) then we would obtain the corresponding additional 
families of densities. 
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2.3 Generalized gamma family emanating 
from equation (1) 

When a — > 1 the forms in (1) for a < 1 and for a > 1 reduces to 

fix) = cx^-'e-"^^', x>0. (20) 

This includes generahzed gamma, gamma, exponential, chisquare, WeibuU, 
Maxwell-Boltzmann, Rayleigh, and related models (Mathai, 1993). The 
model in equation (20) is generally known as the generalized gamma model. 
For 7 = 1,5 = 2, — oo <x<oo,a = l, in (20) we have the Gaussian density 
or error curve, 

fi4:{x) — Ci4 exp (— /3x^), — oo < X < oo. (21) 

For 7 = 5, a = 1 in (20) we have the Wcibull density, given by 

fi5{x)^c,,x'-'e-''^\ x>0. (22) 

For 5 = 1, a = 1 in (20) we have the gamma density, given by 

/i6(a;) = Cl6a:^-V^^ x > 0. (23) 

For 7=1 and /3 — ^ in /ig we have the chisquare density with n degrees of 
freedom, given by 

f^^(x) = cirxt-ig-'/^ a: > 0. (24) 

For 7 = n, /3 = > 0, n = 1,2, ...,6 = 2,a = 1 in (20) we have the chi 

density, given by 

fis{x) = cisx''-'e-^, x>0. (25) 
For 7 = 2, ... in /i6 we have the Erlang density, given by 

/i9 = Cl9xf-le-'^^p = 2, . . . , X > 0. (26) 

For p = 1 in fig we have the exponential density, given by 

f2o{x) = C2oe-^^ a; > 0. (27) 

For 5 = 2, 7 = 3, a = 1 in (20) we have the Maxwell-Boltzmann density, 
given by 

f2i{x) = C2ix^e-^''\ a; > 0. (28) 
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For 5 = 2,7 = 2,/9=2^,o'>0,a = lin (20) we have the Rayleigh density 
(Tirnakh et al., 1998), given by 

f22ix) = C22Xe ^,x>0. (29) 

Remark 2.4 Observe that in (20) if we replace 5 > by —6, 6 > still 
the family belongs to the generalized gamma family. Hence a sequence of 
densities are available by replacing x by - in (20). 

Remark 2.5 If x is replaced by |a;| for —oo < a; < cxo in the models (1), (12) 
and (20) we obtain a series of other densities which will also include double 
exponential or the Laplace density, among others. 

Remark 2.6 In many practical problems there may be a location parameter 
or the variable may be located at a point different from zero. All such cases 
can be covered by replacing x in the pathway model in (1) by (x — b) when 
X — b > or hj \{x — b)\ when — oo < x < oo, where 6 is a constant, called 
the location parameter. 

Observe that in (12) and (20), ^ also belongs to the same family of den- 
sities and hence in (12) and (20) one could have also taken x~^ with 6 > 0. 



3 Pathway model from a generalized entropy 
measure 

A generalized entropy measure of order a, explored by Mathai and Haubold 
(2007), is a generalization of Shannon entropy and it is a variant of the 
generalized entropy of order a in Mathai and Rathie (1975). In the discrete 
case the measure is the following: Consider a multinomial population P — 
(pi, . . . ,Pk), Pi > 0, i = 1, . . . , k, pi + . . . + pk = I. Define the function 

Mk,a{P) = \ « 7^ 1, -oo < a < 2 (30) 

a — 1 

k 

l]m Mk,a{P) = -EPilnp, = Sk{P) (31) 

a— »1 . — ; 



7 



by using L'Hospital's rule. In this notation OlnO is taken as zero when any 
Pi = 0. Thus (30) is a generahzation of Shannon entropy Sk{P) as seen from 
(31). Note that (30) is a variant of Harvda-Charvat entropy H^^aiP) and 
TsaUis entropy Tk,a{P) where 

= a^l,a>0 (32) 

and 

n,a{P) = ^^^^^^, a^l, a>0. (33) 
1 — a 

We will introduce another measure associated with (31) and parallel to Renyi 
entropy Rk,a in the following form: 

In (Ef-ipM 

M*,JP)^ \ a^l,-oo<a<2. (34) 

Renyi entropy is given by 

In ( Y.^- P") 

RkAP) = \ a ^ 1, a > 0. (35) 

1 — a 

It will be seen later that the form in (30) is amenable to the pathway model. 



3.1 Continuous analogue 

The continuous analogue to the measure in (30) is the following: 

MM) ^ ^-'^'-"7-^ (36) 

a — 1 

i-jmr'-mdx - 1 _ E[f{x)r- - 1 

a — 1 a — 1 

where E[-] denotes the expected value of [•]. Note that when a = 1, 

It is easy to see that the generalized entropy measure in (30) is connected 
to Kerridge's "inaccuracy" measure (Kerridge, 1961). The generalized in- 
accuracy measure is E[q{x)Y~°' where the experimenter has assigned q{x) 
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for the true density f{x), where q{x) could be an estimate of f{x) or q{x) 
could be coming from observations. Through disturbance or distortion if the 
experimenter assigns [f{x)]^~°' for then the inaccuracy measure is 

M,(/) of (36). 

3.2 Distributions with maximum generalized entropy 

Among all densities, which one will give a maximum value for Ma{f) in 
equation (36)? Consider all possible functions f{x) such that f{x) > for 
all X, f{x) = outside (a, 6), a < b, /(a) is the same for all such f{x), f{h) 
is the same for all such /, f(x)dx < oo. Let f{x) be a continuous function 
of x with continuous derivatives in (a, b). Let us maximize J^[f{x)Y~°'dx for 
fixed a and over all functional /, under the conditions that the following two 
moment-like expressions be fixed quantities: 

t x^^-^^^^-"^ f{x)dx = given, and f x^'^''^^^^-'^^^^ f{x)dx = given (37) 

Ja J a 

for fixed 7 > and 5 > 0. Consider 

U = [/(x)]'-"-Aix(^-^)(^-'^)/(^) + A2x(^-')^'-")+'^/(x), a < 2, a 7^ 1 (38) 

where Ai and A2 are Lagrangian multipliers. Then the Euler equation is the 
following: 

1^ = ^ (2-a)[/(a;)]i-"-Aix(^-i)M+A2x(^-^)(^-")+'5=0 

^ [/(x)]^-" = -^xM(^-")[l-^a;^] (39) 
(z — a) Ai 

f(x) = c x^-^ [1 - 77(1 - a)x^] ^ (40) 

where X1/X2 is written as 77(1—0;) with 77 > such that 1—77(1— a)a;'^ > since 
f{x) is assumed to be non-negative. By using the conditions (37) and (39) 
wc can determine c and 77. When the range of x for which f{x) is nonzero is 
(0, 00) and when c is a normalizing constant, then (40) is the pathway model 
of Mathai (2005) in the scalar case where a is the pathway parameter. When 
7 = 1, (5 = 1 in (40) then (40) produces the power law. The form in (39) for 
various values of Ai and A2 can produce all the four forms 

aix^-^[l - /?i(l - a)x^]-^, a2X^-^[l - /32(1 - a)x^]^iora < 1 



9 



and 



a3X^-^[l + (3^{a - l)x^\~^, a^x'^-^l + (3i{a - l)x'^]^forQ; > 1 

with ttj, /3j > 0, i = 1, 2, 3, 4. But out of these, the second and the third forms 
can produce densities in (0, oo). The first and fourth will not be converging. 
When f{x) is a density in (40), what is the normalizing constant c? We 
need to consider three cases of a < 1,q; > 1 and a ^ 1. This c is already 
evaluated in section 1. 

Remark 3.1 In Mathai and Haubold (2007) further results are provided 
for the pathway model associated with TsaUis' entropy, fractional calculus, 
Mittag-Leffler functions, and distributions with one of the parameters having 
a prior distribution of its own, giving rise to superstatistics (Beck and Cohen, 
2003; Beck, 2006). 



4 Conclusions 

Based on the pathway model developed by Mathai (2005) we derived a dis- 
tributional pathway proceeding from the generalized typc-1 beta family to 
generalized typc-2 beta family in eq. (12) to generalized gamma family in 
eq. (20). This distributional pathway encompasses, among others, the distri- 
butions of Maxwell- Boltzmann and TsaUis (q-exponential function, see Gell- 
Mann and TsaUis 2004; TsaUis, 2004) that are fundamental to Boltzmann- 
Gibbs statistical mechanics and its nonextensive generalization by TsaUis, 
respectively, as well as Bose-Einstein and Fermi-Dirac distributions. Subse- 
quent to the distributional pathway, an entropic pathway, emanates from eq. 
(1) by maximum principle applied to the generalized entropic form of order a 
in (30), covering entropic functionals of Shannon, Boltzmann-Gibbs, Renyi, 
TsaUis, and Harvda-Charvat (Mathai and Haubold, 2007). The results in 
this paper are a contribution to the on-going debate in the physical litera- 
ture concerning the generalization of Boltzmann-Gibbs statistical mechanics 
through a one-parameter generalization of Boltzmann-Gibbs entropy, known 
as TsaUis entropy contained in the seminal paper of TsaUis (1988), whose 
form is equivalent to that of Harvda-Charvat, and which opened the door to 
generalizations of Boltzmann-Gibbs statistical mechanics. 
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